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Abstract 

We address electron transport in honeycomb lattice ribbons with armchair edges attached to two semi- 
infinite one-dimensional metallic electrodes within the tight-binding framework. Here we present numer- 
ically the conductance-energy and current-voltage characteristics as functions of the length and width 
of the ribbons. Our theoretical results predict that for a ribbon with much smaller length and width, 
so-called a nanoribbon, a gap in the conductance spectrum appears across the energy E = 0. While, this 
gap decreases gradually with the increase of the size of the ribbon, and eventually it almost vanishes. 
This reveals a transformation from the semiconducting to the conducting material, and it becomes much 
more clearly visible from our presented current-voltage characteristics. 
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1 Introduction 

The electronic transport in nanoribbons of 
graphcne has opened up new areas in nanoelectron- 
ics. A graphene nanoribbon (GNR) is a monolayer 
of carbon atoms arranged in a honeycomb lattice 
structure [JJ [3J 03 H] . Due to the special electronic 
and physical properties, graphene based materials 
exhibit several novel properties like unconventional 
quantum Hall effect [5], high carrier mobility [3] 
and many others. The high carrier mobility in 
graphene demonstrates the idea for fabrication of 
high speed switching devices those have widespread 
applications in different fields. Some recent experi- 
ments [6j [8] have also suggested that GNRs can 
be used to design field-effect transistors and this 
application provides a huge interest in the commu- 
nity of nanoelectronics device research. Further- 
more, GNRs can be used to construct MOSFETs 
which perform much better than conventional Si 
MOSFETs. In other experiment ,9] it has been pro- 
posed that a narrow strip of graphene with arm- 
chair edges, so-called a graphene nanoribbon, ex- 
hibits semiconducting behavior due to its edge ef- 
fects, unlike carbon nanotubes of larger sizes which 
are mixtures of both metallic and semiconducting 
materials. This is due to the fact that in a nar- 
row graphene sheet, a band gap appears across the 
energy E — 0, while the gap gradually disappears 
with the increase of the size of the ribbon. This 
reveals a transformation from the semiconducting 
to the metallic material, and such a phenomenon 
can be utilized for fabrication of electronic devices. 
This motivates us to study the electron transport in 
honeycomb lattice ribbons with armchair edges and 
to verify qualitatively how the transformation from 
the semiconducting to the conducting property can 
be achieved simply by tuning the size of a ribbon. 

The purpose of the present paper is to provide 
a qualitative study of electron transport in honey- 
comb lattice ribbons with armchair edges attached 
to two semi-infinite one-dimensional metallic elec- 
trodes (see Fig. [TJ). The theoretical description 
of electron transport in a bridge system has been 
followed based on the pioneering work of Aviram 
and Ratner [lOj . Later, many excellent experi- 
ments [TTJ H21 H3 HH US] have been done in several 
bridge systems to understand the basic mechanisms 
underlying the electron transport. Though in liter- 
ature many theoretical [TBI EZl [TBI HH1 HOI I2H 123 H3 
[2l[2g[2S[2a[2l[2S[^[3ll[32]as well as experi- 
mental papers [JTJ [12j [13l [Til HI] on electron trans- 
port are available, yet lot of controversies are still 



present between the theory and experiment, and the 
complete knowledge of the conduction mechanism 
in this scale is not very well established even to- 
day. Several controlling parameters are there which 
can regulate significantly the electron transport in 
a conducting bridge, and all these effects have to be 
taken into account properly to reveal the transport 
properties. For our illustrative purposes, here we 
describe very briefly some of these effects. 

(i) The quantum interference effect [33[ [34] [35j [36j 
13 T] of electron waves passing through different arms 
of any conducting element which bridges two elec- 
trodes becomes the most significant issue. 

(ii) The coupling of the electrodes with bridging 
material provides an important signature in the de- 
termination of current amplitude across any bridge 
system [33]. The understanding of this coupling to 
the electrodes under non-equilibrium condition is a 
major challenge, and we should take care about it 
in fabrication of any electronic device. 

(iii) The geometry of the conducting material be- 
tween the two electrodes itself is an important issue 
to control the electron transmission. To emphasize 
it, Ernzerhof et al. 38J have predicted several model 
calculations and provided some significant results. 

(iv) The dynamical fluctuation in the small-scale 
devices is another important factor which plays 
an active role and can be manifested through the 
measurement of shot noise [39[ 140] . a direct con- 
sequence of the quantization of charge. It can be 
used to obtain information on a system which is 
not available directly through the conductance mea- 
surements, and is generally more sensitive to the 
effects of electron-electron correlations than the av- 
erage conductance. 

Furthermore, several other parameters of the 
Hamiltonian that describe a system also provide sig- 
nificant effects in the determination of the current 
across a bridge system. 

Here we adopt a simple tight-binding model to 
describe the system and all the calculations are per- 
formed numerically. We address the conductance- 
energy and current-voltage characteristics as func- 
tions of lengths and widths of ribbons. Our results 
clearly predicts how a honeycomb lattice ribbon 
with armchair edges transforms its behavior from 
the semiconducting to the metallic nature, and this 
feature may be utilized in fabrication of nanoelec- 
tronic devices. 

The paper is organized as follow. Following the 
introduction (Section 1), in Section 2, we present 
the model and the theoretical formulations for our 
calculations. Section 3 discusses the significant re- 
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suits, and finally, we summarize our results in Sec- full system i.e., the ribbon, source and drain, the 
tion 4. Green's function is defined as, 



2 Model and the synopsis of 
the theoretical background 

Let us refer to Fig. Q] where a honeycomb lat- 
tice ribbon with armchair edges is attached to two 
semi-infinite one-dimensional metallic electrodes, 
viz, source and drain. It is important to note that 
throughout this study we attach the electrodes at 




Honeycomb lattice ribbon 



Drain 



Figure 1: Schematic view of a honeycomb lattice 
ribbon with armchair edges attached to two semi- 
infinite one-dimensional metallic electrodes, viz, 
source and drain. Filled circles correspond to the 
position of the atomic sites (for color illustration, 
see the web version). 

the two extreme ends of nanoribbons, as seen in 
Fig. [TJ to keep the uniformity of the quantum in- 
terference effects. 

To calculate the conductance g of the ribbon, we 
use the Landauer conductance formula 41, 42], and 
at very low temperature and bias voltage it can be 
expressed in the form, 



2e 2 



-T 



(1) 



where T gives the transmission probability of an 
electron in the ribbon. This (T) can be represented 
in terms of the Green's function of the ribbon and 
its coupling to the two electrodes by the relation [41j 



T = Tr \T s G r rib T D G a rib ] 



(2) 



where G r rib and G° ib are respectively the retarded 
and advanced Green's functions of the ribbon in- 
cluding the effects of the electrodes. The parame- 
ters V $ and Tjj describe the coupling of the ribbon 
to the source and drain respectively, and they can 
be defined in terms of their self-energies. For the 



G=(e-HX 



(3) 



where e = E + i5. E is the injecting energy of the 
source electron and 5 gives an infinitesimal imagi- 
nary part to e. To Evaluate this Green's function, 
the inversion of an infinite matrix is needed since 
the full system consists of the finite ribbon and the 
two semi-infinite electrodes. However, the entire 
system can be partitioned into sub-matrices cor- 
responding to the individual sub-systems and the 
Green's function for the ribbon can be effectively 
written as, 



G r ib — (e — H, 



rib 



TiQ — E; 



(4) 



where H r i b is the Hamiltonian of the ribbon which 
can be written in the tight-binding model within 
the non-interacting picture like, 



Hrib — ^ y C i Ci 



E 

<ij> 



(4 



(5) 



In the above Hamiltonian (H r i b ), e^'s are the site 
energies, cj (a) is the creation (annihilation) oper- 
ator of an electron at the site i and t is the nearest- 
neighbor hopping integral. Similar kind of tight- 
binding Hamiltonian is also used to describe the 
two semi-infinite one-dimensional perfect electrodes 
where the Hamiltonian is parametrized by constant 
on-site potential eo and nearest-neighbor hopping 
integral t . In Eq. 01 E s = h s _ rib g s h S - r ib and 
= h D - rib g D W D _ rib are the self-energy opera- 
tors due to the two electrodes, where gs and gn 
correspond to the Green's functions of the source 
and drain respectively. hs- r ib and ho-rib are the 
coupling matrices and they will be non-zero only 
for the adjacent points of the ribbon, and the elec- 
trodes respectively. The matrices Ts and Tjj can 
be calculated through the expression, 



r 



S(D) 



J S(D) 



2j S(D) 



(6) 



where and are the retarded and ad- 

vanced self-energies respectively, and they are con- 
jugate with each other. These self-energies can be 
written as [43] . 



J S{D) 



A 



S(D) 



iA 



S(D) 



(7) 



where are the real parts of the self-energies 

which correspond to the shift of the energy eigenval- 
ues of the ribbon and the imaginary parts A s(£>) of 
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the self-energies represent the broadening of these 
energy levels. Since this broadening is much larger 
than the thermal broadening, we restrict our all cal- 
culations only at absolute zero temperature. All 
the informations about the ribbon-to-electrode cou- 
pling arc included into these two self-energies. 

The current passing across the ribbon can be de- 
picted as a single-electron scattering process be- 
tween the two reservoirs of charge carriers. The 
current I can be computed as a function of the ap- 
plied bias voltage V through the relation [41], 



I(V) 



e 

■kTi 



E F +eV/2 



E F -eV/2 



T(E, V)dE 



(8) 



where Ep is the equilibrium Fermi energy. Here we 
make a realistic assumption that the entire voltage 
is dropped across the ribbon-electrode interfaces, 
and it is examined that under such an assumption 
the I-V characteristics do not change their qualita- 
tive features. This assumption is based on the fact 
that, the electric field inside the ribbon especially 
for narrow ribbons seems to have a minimal effect 
on the conductance-voltage characteristics. On the 
other hand, for quite larger ribbons and high bias 
voltages the electric field inside the ribbon may play 
a more significant role depending on the internal 
structure and size of the ribbon [43] , but the effect 
becomes too small. 

3 Results and discussion 

In order to understand the dependence of electron 
transport on the lengths and widths of nanorib- 
bons, in the present article, we concentrate only 
on the cleaned systems rather than any dirty one. 
Accordingly, we set the site energies of the honey- 
comb lattice ribbons as = for all i. The values 
of the other parameters are assigned as follow: the 
nearest-neighbor hopping integral t in the ribbon is 
set to 2, the on-site energy eo and the hopping in- 
tegral to for the two electrodes are fixed to and 2 
respectively. The parameters t$ and Tp> are set as 
1.5, where they correspond to the hopping strengths 
of the ribbon to the source and drain respectively. 
In addition to these, we also introduce two other 
parameters N and M to reveal the size of a nanorib- 
bon, where they correspond to the width and length 
of the ribbon respectively. Thus, for example, a 
nanoribbon with N = 1 and M = 4 represents a 
linear chain of four hexagons. Hence the parame- 
ter M determines the total number of hexagons in 
a single chain. Following this rule, a nanoribbon 



with N = 3 and M = 3 corresponds to three lin- 
ear chains attached side by side (see Fig. []} where 
each chain contains three hexagons. For simplic- 
ity, throughout our study we set the Fermi energy 
Ep = and choose the units where c = e = h = 1. 

Let us first describe the variation of the conduc- 
tance g as a function of the injecting electron energy 
E. In Fig. Owe present the conductance-energy (g- 
E) characteristics for some honeycomb lattice rib- 
bons with fixed width (iV = 1) and varying lengths, 
where (a) and (b) correspond to the linear chains 
with six (M = 6) and ten (M = 10) hexagons 
respectively. The conductance spectra shows fine 
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Figure 2: Conductance g as a function of the energy 
E for some lattice ribbons with fixed width N = 1 
and varying lengths where (a) M = 6 and (b) M — 
10 (for color illustration, see the web version). 

resonance peaks for some particular energies, while 
for all other values of the energy E, either it (g) 
drops to zero or gets much small value. At these res- 
onance energies, the conductance gets the value 2, 
and hence, the transmission probability T becomes 
unity since the expression g — 2T holds from the 
Landauer conductance formula (see Eq. [T]) . These 
resonance peaks are associated with the energy lev- 
els of the nanoribbons and thus the conductance 
spectra, on the other hand, reveal the signature of 
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the energy spectra of the nanoribbons. The most 
important issue observed from these spectra is that, 
a central gap appears across the energy E — and 
the width of the gap becomes small for the chain 
with 10 hexagons compared to the other chain i.e., 
the chain with 6 hexagons. It predicts that, for a 
fixed width, the central energy gap decreases with 
the increase of the length of the nanoribbon. In 
the same footing, to visualize the dependence of the 
width on the conductance-energy characteristics, in 
Fig. [3] we display the results for some honeycomb 
lattice ribbons considering the width TV = 4, where 




Figure 3: Conductance g as a function of the energy 
E for some lattice ribbons with fixed width N = 4 
and varying lengths (identical as in Fig. [2]) where 
(a) M = 6 and (b) M = 10. Here the width of 
the ribbons is increased compared to the ribbons as 
taken in Fig. [2] (for color illustration, see the web 
version) . 

(a) and (b) represent the nanoribbons with identi- 
cal lengths as in Fig. [2l The results show that, due 
to the large system sizes the g-E characteristics ex- 
hibit almost a quasi-continuous variation across the 
energy E — 0. For both these two ribbons the en- 
ergy gap also appears around the energy E = 0, and 
the gap decreases with the increase of the length of 
the nanoribbon. Comparing the results presented 



in Figs. [2] and [3l we can emphasize that for a fixed 
width the central energy gap always decreases with 
the size of the nanoribbon. Now to reveal the de- 
pendence of the energy gap on the system size much 
more clearly, in Fig. [4] we show the variation of the 
central energy gap 6E as a function of the length 
M for some honeycomb lattice ribbons with differ- 
ent widths N. The red, green and blue lines cor- 
respond to the results for the ribbons with fixed 
widths N = 1, 2 and 4 respectively. These results 
clearly emphasize that for the fixed width the gap 
gradually decreases with the increase of the length 
of the nanoribbon. It is also examined that for much 
larger lengths it (SE) almost vanishes (not shown 
here in the figure). Quite similar nature is also ob- 
served if we plot the variation of the energy gap as 
a function of the length N keeping the width M 
as a constant, and due to the obvious reason we 
do not plot these results further in the present de- 
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Figure 4: Variation of the central energy gap SE as 
a function of the length M for some lattice ribbons 
with fixed widths N. The red, green and blue curves 
correspond to N = 1, 2 and 4 respectively (for color 
illustration, see the web version). 

scription. These results provide us an important 
signature which concern with a transition from the 
semiconducting (finite energy gap) to the conduct- 
ing (zero energy gap) material, and this transition 
can be achieved simply by tuning the size of the 
nanoribbon. 

All these basic features of electron transfer can 
be quite easily explained from our study of the 
current-voltage (I-V) characteristics rather than 
the conductance-energy spectra. The current / is 
determined from the integration procedure of the 
transmission function (T) (see Eq. [8J, where the 
function T varies exactly similar to the conduc- 
tance spectra, differ only in magnitude by a fac- 
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tor 2, since the relation g = 2T holds from the 
Landauer conductance formula (Eq. [1]) . The varia- 
tion of the current-voltage characteristics for some 
typical honeycomb lattice ribbons with fixed width 
N = 2 and varying lengths is presented in Fig. 
where (a) and (b) correspond to the ribbons with 
M = 3 and 5 respectively. The current exhibits 
a staircase like behavior as a function of the ap- 



strong coupling limit, described by the condition 
T S(D) ~ t, current varies quite continuously with 
the bias voltage V and achieves large current am- 
plitude compared to the weak-coupling limit. All 
these coupling effects have clearly been explained 
in many papers in the literature. The significant 
feature observed from the figure (Fig. [5|) is that for 
the fixed width (TV = 2), the threshold bias volt- 
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Figure 5: Current I as a function of the bias voltage 
V for some lattice ribbons with fixed width N = 2 
and varying lengths where (a) M = 3 and (b) M — 
5 (for color illustration, see the web version). 

plied bias voltage V. This staircase like nature ap- 
pears due to the existence of the resonance peaks in 
the conductance spectra since the current is com- 
puted by the integration process of the transmis- 
sion function T. As we increase the bias voltage 
V, the electrochemical potentials in the two elec- 
trodes cross one of the energy levels of the ribbon 
and accordingly a jump in the I-V curve appears. 
The sharpness of the steps in the current-voltage 
characteristics and the current amplitude solely de- 
pend on the coupling strengths of the nanoribbon 
to the electrodes, viz, source and drain. It is ob- 
served that, in the limit of weak coupling, defined 
by the condition t$(d) « t, current shows stair- 
case like structure with sharp steps. While, in the 



Figure 6: Current I as a function of the bias voltage 
V for some lattice ribbons with fixed width N = 3 
and varying lengths where (a) M = 2 and (b) M = 
3 (for color illustration, see the web version). 

age (Vth) of electron conduction decreases with the 
increase of the length of the ribbon. This reveals 
a transformation towards the conducting material. 
Quite in the same fashion, to see the variation of the 
threshold bias voltage Vth for other system sizes, in 
Fig. [5] we plot the results for some nanoribbons with 
fixed width N = 3 and varying lengths where (a) 
and (b) correspond to the ribbons with M — 2 and 
3 respectively. The results show that the thresh- 
old bias voltages decrease much more compared to 
the nanoribbons of width N = 2. Thus both from 
Figs. and [S] we clearly observe that Vth can be reg- 
ulated very nicely by tuning the size (both length 
and width) of the nanoribbon. For quite larger rib- 
bons the threshold bias voltage eventually reduces 
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to zero. This clearly manifests the transformation 
from the semiconducting to the conducting mate- 
rial. 

4 Concluding remarks 

To summarize, we have addressed electron trans- 
port in honeycomb lattice ribbons with armchair 
edges attached to two semi-infinite one-dimensional 
metallic electrodes within the tight-binding frame- 
work. Our results have predicted that for ribbons 
with smaller lengths and widths, a central gap in 
the conductance spectrum appears across the en- 
ergy E = 0. While, this gap decreases gradually as 
we increase the size of the ribbon and eventually it 
almost vanishes. This reveals a transformation from 
the semiconducting to the conducting behavior, and 
it has been much more clearly described from the 
presented current-voltage characteristics which pro- 
vide that the threshold bias voltage of electron con- 
duction decreases gradually with the increase of the 
size of the ribbon. 

This is our first step to describe how the elec- 
tron transport in honeycomb lattice ribbons with 
armchair edges depends on the length and width of 
the ribbons. Here we have made several realistic 
assumptions by ignoring the effects of the electron- 
electron correlation, disorder, temperature, etc. All 
these effects can be incorporated quite easily with 
this present formalism and we need further study in 
such systems. In our next work we are investigat- 
ing the electron transport properties in honeycomb 
lattice ribbons with zigzag edges. 
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